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Abstract. We study the structure of newly born neutron stars (protoneutron stars) within the finite temperature Brueckner- 
Bethe-Goldstone theoretical approach including also hyperons. We find that for purely nucleonic stars both finite temperature 
and neutrino trapping reduce the value of the maximum mass. For hyperonic stars the effect is reversed, because neutrino 
trapping shifts the appearance of hyperons to larger baryon density and stiffens considerably the equation of state. 
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^ 1. Introduction 

> : 

A protoneutron star (PNS) is formed after a successful su- 
I pernova explosion and constitutes for several tens of seconds 
a transitional state to either a neutron star or a black hole 
<3 (Prakash et al. 1997). Initially, the PNS is optically thick to neu- 
trinos, that is, they are temporarily trapped within the star. The 
| subsequent evolution of the PNS is dominated by neutrino dif- 
r"| fusion, which first results in deleptonization and subsequently 
"^T* | m cooling. After a much longer time, photon emission com- 
f q ■ petes with neutrino emission in neutron star cooling. 

In this paper, we will focus upon the essential ingredient 
Ph ■ that governs the macrophysical evolution of neutron stars, i.e., 
;> the equation of state (EOS) of dense matter at finite temper- 
k> ' ature - We have developed a microscopic EOS in the frame- 
r% \ work of the Brueckner-Bethe-Goldstone (BBG) many-body ap- 
■ proach extended to finite temperature. This EOS has been suc- 
cessfully applied to the study of the limiting temperature in 
nuclei (Baldo et al. 1999, 2004). The scope of this work is to 
present results on the composition and structure of these newly 
born stars with the EOS previously mentioned. 

Two new effects have to be considered for a PNS. First, 
thermal effects which result in entropy production, with values 
of few units per baryon, and temperatures up to 30-40 MeV 
(Burrows & Lattimer 1986; Pons et al. 1999). Second, the fact 
that neutrinos are trapped in the star, which means that the neu- 
trino chemical potential is non-zero. This alters the chemical 
equilibrium and leads to compositional changes. Both effects 
may result in observable consequences in the neutrino signa- 
ture from a supernova and may also play an important role in 
determining whether or not a given supernova ultimately pro- 
duces a cold neutron star or a black hole. 

This paper is organized as follows. In Section 2 we briefly 
illustrate the BBG many-body theory at finite temperature. 
Section 3 is devoted to the study of the stellar matter composi- 



tion, both with and without neutrino trapping, and the resulting 
EOS. In Section 4 we discuss our results regarding the struc- 
ture of (proto)neutron stars, in particular their maximum mass. 
Finally, conclusions are drawn in Section 5. 

2. Brueckner theory 

Over the last two decades the increasing interest for the EOS 
of nuclear matter has stimulated a great deal of theoretical ac- 
tivity. Phenomenological and microscopic models of the EOS 
have been developed along parallel lines with complementary 
roles. The former models include nonrelativistic mean field the- 
ory based on Skyrme interactions (Bone he et al. 1998) and rel- 
ativistic mean field theory based on meson-exchange interac- 
tions (Walecka model) (Serot & Walecka 1986). Both of them 
fit the parameters of the interaction in order to reproduce the 
empirical saturation properties of nuclear matter extracted from 
the nuclear mass table. The latter ones include nonrelativis- 
tic Brueckner-Hartree-Fock (BHF) theory (Baldo 1999) and its 
relativistic counterpart, the Dirac-Brueckner theory (Machleidt 
1989; Li et al. 1992), the nonrelativistic variational approach 
also corrected by relativistic effects (Akmal et al. 1997, 1998; 
Morales et al. 2002), and more recently the chiral perturbation 
theory (Kaiser et al. 2002). In these approaches the parame- 
ters of the interaction are fixed by the experimental nucleon- 
nucleon and/or nucleon-meson scattering data. 

For states of nuclear matter with high density and high 
isospin asymmetry the experimental constraints on the EOS are 
rather scarse and indirect. Different approaches lead to differ- 
ent or even contradictory theoretical predictions for the nuclear 
matter properties. The interest for these properties lies, to a 
large extent, in the study of dense astrophysical objects, i.e., 
supernovae and neutron stars. 

One of the most advanced microscopic approaches to the 
EOS of nuclear matter is the Brueckner theory. In the recent 



2 O. E. Nicotra, M. Baldo, G. F. Burgio, and H.-J. Schulze: Protoneutron stars within the Brueckner-Bethe-Goldstone theory 



years, it has made a rapid progress in several aspects: (i) The 
convergence of the Brueckner-Bethe-Goldstone expansion has 
been firmly established (Song et al. 1998; Baldo et al. 2000b). 
(ii) The addition of phenomenological three-body forces (TBF) 
based on the Urbana model (Carlson et al. 1983; Schiavilla et 
al. 1986), permitted to improve to a large extent the agreement 
with the empirical saturation properties (Baldo et al. 1997; 
Zhou et al. 2004). (iii) The BHF approach has been extended 
in a fully microscopic and self-consistent way to describe nu- 
clear matter containing also hyperons (Schulze et al. 1995, 

1998) , opening new fields of applications such as hypernuclei 
(Cugnon et al. 2000; Vidana et al. 2001) and a more realistic 
modeling of neutron stars (Baldo et al. 1998, 2000a). 

In this paper, we will discuss the nuclear EOS at finite tem- 
perature within the BBG theory, and apply it to the study of 
protoneutron stars. For this purpose, we have performed a sys- 
tematic calculation of the EOS within the Bloch-De Dominicis 
formalism with a realistic nucleon-nucleon interaction, includ- 
ing nucleonic three-body forces. Details are given in the fol- 
lowing section. 

2.1. EOS of nuclear matter at finite temperature 

Only few microscopic calculations of the nuclear EOS at fi- 
nite temperature are so far available. The variational calcu- 
lation by Friedman & Pandharipande (1981) was one of the 
first semi-microscopic investigations of the finite temperature 
EOS. The results predict a Van der Waals behavior, which 
leads to a liquid-gas phase transition with a critical temperature 
T c =s 18-20 MeV. Later, Brueckner calculations (Lejeune et al. 
1986; Baldo & Ferreira 1999) and chiral perturbation theory at 
finite temperature (Kaiser et al. 2002) confirmed these findings 
with very similar values of T c . The Van der Waals behavior was 
also found in the finite temperature relativistic Dirac -Brueckner 
calculations of (Ter Haar & Malfliet 1986, 1987; Huber et al. 

1999) , although at a lower temperature. 

The formalism which is closest to the BBG expansion, and 
actually reduces to it in the zero temperature limit, is the one 
formulated by Bloch & De Dominicis (1958, 1959a, 1959b). In 
this approach the essential ingredient is the two-body scatter- 
ing matrix K, which, along with the single-particle potential U, 
satisfies the self-consistent equations 



{k x k 2 \K{W)\hk A ) = (hk 2 \V\hh) 



+ Re Yihhmk'k') 



[l-n(k' 3 )][l-n(k\)] 
W-E K - E K + ie 



{k' 3 k' 4 \K(W)\k 3 k 4 ) 



(1) 



and 



U(k 1 )=Yjn(k2)(k 1 k 2 \K(W)\k 1 k 2 ) A , 



(2) 



where kj generally denote momentum, spin, and isospin. Here 
V is the two-body interaction, W = + E^ represents the 
starting energy, and Et = k 2 /2m + U(k) the single-particle 
energy. Eq. Q coincides with the Brueckner equation for the 
K matrix at zero temperature, if the single-particle occupation 
numbers n(k) are taken at T = 0. At finite temperature n(k) 



is a Fermi distribution. For a given density and temperature, 
Eqs. Q and (0 have to be solved self-consistently along with 
the following equation for the auxiliary chemical potential ft, 



(3) 



The grand-canonical potential density to can be written (Baldo 
& Ferreira 1999), 

u> = c0{) - n(k)U(k) 

k 

+ l - J] n(h)n{k 2 ) {hk 2 \K{W)\k,k 2 ) A , (4) 

kite 

where a>o is the grand-canonical potential density for a gas of 
independent particles with single-particle spectrum Ej, 



0)Q 



-IVln(l 



) 



(5) 



Eq. 10} neglects a series of terms proportional to n(k)[l - n(k)] 
(or powers of it), which turn out to be negligible in the temper- 
ature and density ranges relevant for neutron and protoneutron 
stars (Baldo 1999). The free energy density is then 



/ = (o + pp. 



(6) 



and the "true" chemical potential p and the pressure p are given 
by 



df 
dp 



P = P 



, d(f/p) 

dp 



(7) 
(8) 



We stress that this procedure permits to fulfill the Hugenholtz- 
Van Hove theorem in the calculation of thermodynamic al quan- 
tities in the Brueckner theory. For an extensive discussion of 
this topic, the reader is referred to (Baldo 1999), and references 
therein. 

The determination of the different configurations inside 
protoneutron stars requires the calculation of the EOS at dif- 
ferent chemical compositions in search of the beta equilibrium. 
To save computational time and simplify the numerical work 
we introduce an approximate procedure. For each configuration 
the single-particle potential Ui(k) for the component i is calcu- 
lated self-consistently at T — 0. At T + we consider £/,(£) in- 
dependent of T. In this approximation the correlations at T + 
are assumed to be essentially the same as at T = {Frozen 
Correlations Approximation). It turns out that the assumed in- 
dependence is valid to a good accuracy (Baldo & Ferreira 1999, 
Fig. 12), at least for not too high temperature, due to a substan- 
tial compensation between the T-dependence of the A'-matrix 
and the T-dependence of the Fermi distributions, see Eq. (J3- 
Once this approximation is introduced, the grand-canonical po- 
tential density to can be trivially calculated and the free energy 
density has the following simplified expression 



(9) 
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where 

si=~Yi ( ni{k) ]nni(k) + [1 - ni(k)] ln[1 - ni(k)] ) (10) 

k 

is the entropy density for component i treated as a free gas with 
spectrum Ej(k). To illustrate the accuracy of this approxima- 
tion, in Fig.^is reported the comparison between the Frozen 
Correlations Approximation and the full microscopic calcula- 
tions for the free energy in symmetric nuclear matter in the 
relevant density range. For pure neutron matter the agreement 
is even closer. 



2.2. Effects of three-body forces 

It is well known that at zero temperature the non-relativistic 
microscopic approaches do not reproduce correctly the nuclear 
matter saturation point po a 0.17 fnT 3 , E/A a -16 MeV. One 
common way of correcting this deficiency is to introduce three- 
body forces among nucleons. A complete microscopic theory 
of TBF is not available yet, thus one is forced to use phe- 
nomenological models. We have adopted the so-called Urbana 
model (Carlson et al. 1983; Schiavilla et al. 1986), which con- 
sists of an attractive term due to two-pion exchange with ex- 
citation of an intermediate A resonance, and a repulsive phe- 
nomenological central term. In the BBG approach the TBF 



is reduced to a density dependent two-body force by averag- 
ing over the position of the third particle, assuming that the 
probability of having two particles at a given distance is re- 
duced according to the two-body correlation function. The cor- 
responding EOS satisfies several requirements, namely, (i) it 
reproduces correctly the nuclear matter saturation point (Baldo 
et al. 1997; Zhou et al. 2004), (ii) the incompressibility is com- 
patible with values extracted from phenomenology (Myers & 
Swiatecki 1996), (iii) the symmetry energy is compatible with 
nuclear phenomenology, (iv) the causality condition is fulfilled 
up to the densities typical of neutron star cores. In all calcula- 
tions presented in this paper we use the Argonne Vig nucleon- 
nucleon potential (Wiringa et al. 1995) along with the phe- 
nomenological Urbana TBF. 

In Fig.|2]we display the EOS obtained following the above 
discussed procedure, both for symmetric and purely neutron 
matter. In the upper panels we display the free energy, whereas 
the lower panels show the internal energy per particle as a func- 
tion of the nucleon density, for several values of the tempera- 
ture between and 50 MeV. We notice that the free energy of 
symmetric matter shows a typical Van der Waals behavior (with 
Tc ~ 16 MeV) and is a monotonically decreasing function of 
the temperature. At T — the free energy coincides with the to- 
tal energy and the corresponding curve is just the usual nuclear 
matter saturation curve. On the contrary, the internal energy is 
an increasing function of the temperature. The effect is less pro- 
nounced for pure neutron matter due to the larger Fermi energy 
of the neutrons at given density. 

2.3. Inclusion of hyperons 

It is commonly accepted that, whereas at moderate densities 
p ~ po the matter inside a neutron star consists only of nucle- 
ons and leptons, at higher densities other baryonic species may 
appear due to the fast rise of the baryon chemical potentials 
with density (Glendenning 1982, 1985). For this purpose we 
have extended the BBG approach in order to include the E and 
A hyperons (Schulze et al. 1998; Baldo et al. 1998, 2000a). The 
inclusion of hyperons is a delicate task, and requires the knowl- 
edge of the nucleon-hyperon (NH) and hyperon-hyperon (HH) 
interactions. In our work we have used the Nijmegen soft-core 
NH potential (Maessen et al. 1989), that is well adapted to the 
available experimental NH scattering data. Unfortunately, no 
HH scattering data and therefore no reliable HH potentials are 
available at the moment. We have therefore neglected these in- 
teractions in our calculations, which is supposedly justified, as 
long as the hyperonic partial densities remain limited. 

We have found that due to its negative charge the E~ hy- 
peron is the first strange baryon appearing in the reaction n + 
n — > p+Jr, in spite of its substantially larger mass compared to 
the neutral A hyperon (M z - = 1 197 MeV, M A = 1 1 16 MeV). 
The presence of hyperons strongly softens the EOS, mainly due 
to the larger number of baryonic degrees of freedom. This EOS 
produces a maximum neutron star mass that lies slightly be- 
low the canonical value of 1.44 M G (Taylor & Weisberg 1989), 
which could indicate the presence of non-baryonic (quark) mat- 
ter in the interior of heavy neutron stars (Burgio et al. 2002a, 
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Fig. 2. Finite temperature BBG equation of state for symmetric (left-hand panels) and purely neutron (right-hand panels) matter. 
The upper panels show the free energy and the lower panels the internal energy per particle as a function of the nucleon density. 
The temperatures vary from to 50 MeV in steps of 10 MeV. 



2002b; Baldo et al. 2003, Maieron et al. 2004). However, the 
quantitative effects of more reliable NH and HH potentials as 
well as hyperonic TBF need still to be explored in the future. 

For these reasons, we will present in this article only 
zero-temperature results with interacting hyperons and perform 
finite-temperature calculations in a much simpler way using 
non-interacting hyperons in order to estimate the sensitivity of 
the main global observables to finite temperature. 

3. Composition and EOS of hot stellar matter 

For stars in which the strongly interacting particles are only 
baryons, the composition is determined by the requirements of 
charge neutrality and equilibrium under the weak processes 

Bi -> B 2 + l + v, , B 2 + l^B x +vi, (11) 

where B\ and B2 are baryons and / is a lepton, either an electron 
or a muon. When the neutrinos have left the system, these two 
requirements imply that the relations 

^ qiXj + 2^ qm = (12) 

i 1 



and 

Hi = bifx n - qiHi (13) 

are satisfied. In the expression above, Xj = pi/ps represents the 
baryon fraction for the species i, and pg the baryon density. 
The neutron chemical potential is denoted by //„, whereas Hi 
refers to the chemical potential of the baryon species i, bi to its 
baryon number and qi to the electric charge. The same holds 
true for the quantities with subscript I. Under the condition that 
the neutrinos are trapped in the system, the beta equilibrium 
equation Jl 3t is altered to 

Hi = bifi„ - qiipii - fly) , (14) 

where Hv, is the chemical potential of the neutrino v/. 

For stellar matter containing nucleons and hyperons as rel- 
evant degrees of freedom, the chemical equilibrium conditions 
read explicitly 

H„ - \ip = He- Hv e =Hn+ (*v„ , 
Hz- = 2h„ - Hp . 

Ha = Hn ■ (15) 
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Fig. 3. Relative populations for neutrino-free matter as a function of the baryon density for several values of the temperature. 
Upper panels show the partial concentrations for the purely nucleonic case, whereas in the lower panels noninteracting hyperons 
£~ and A are included as well. 



Because of trapping, the numbers of leptons per baryon of each 
flavor I - e,p, 



(16) 



are conserved on dynamical time scales. Gravitational collapse 
calculations of the white-dwarf core of massive stars indicate 
that at the onset of trapping, the electron lepton number Y e = 
x e + x Ve 0.4, the precise value depending on the efficiency 
of electron capture reactions during the initial collapse stage. 
Also, because no muons are present when neutrinos become 
trapped, the constraint Y M = x M - Xy = can be imposed. We 
fix the Yi at these values in our calculations for neutrino-trapped 
matter. 

The determination of the chemical potentials is one of the 
most important and delicate steps of the entire calculation. In 
fact, in asymmetric nuclear matter, it involves the knowledge 
of the free energy and its partial derivatives with respect to the 
total baryon density and proton fraction, i.e., 



ix n (p,x p ) = 
H P (p,x p ) = 



d d 

1 +P7 X p 



dp 



dx n 



f 



8 ,1 N 9 

1 +p— +(1 'Xp) — 

op ox p \ p 



(17) 
(18) 



The derivatives are computed by using an analytical fit of the 
free energy f(T,p, x p ), which is given in Appendix A. 

As far as the hyperon chemical potentials are concerned, 
they are obtained through the standard procedure at zero tem- 
perature illustrated in (Maieron et al. 2004). The chemical po- 
tentials of the noninteracting leptons and hyperons at finite tem- 
perature are obtained by solving numerically the free Fermi gas 
model at finite temperature. 

Once the baryonic and leptonic chemical potentials are 
known, one can proceed to calculate the composition of the 
/S- stable stellar matter, and then the free energy density / and 
pressure p through the usual thermodynamical relation 

2 d(f/p) 



P=P 



d P 



(19) 



The stable configurations of a (proto)neutron star can be ob- 
tained from the well-known hydrostatic equilibrium equations 
of Tolman, Oppenheimer, and Volkov (Shapiro & Teukolsky 
1983) for the pressure p and the enclosed mass m 

dp(r) Gm(r)e(r) [l + + ^ 



dr 

dm(r) 
dr 



1 



2Cm(r) 



= 4nr 2 e(r) , 



(20) 



(21) 
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Fig. 4. Same as Fig. [3] but for neutrino-trapped matter. 



once the EOS p(e) is specified, being e the total energy density 
(G is the gravitational constant). For a chosen central value of 
the energy density, the numerical integration of Eqs. J20l > and 
( 12 li provides the mass-radius relation. 

For the description of the NS crust at zero temperature, we 
have joined the hadronic equations of state above described 
with the ones by Negele & Vautherin (1973) in the medium- 
density regime (0.001 ftrT 3 < p < 0.08 fm -3 ), and the ones by 
Feynman, Metropolis, & Teller (1949) and Baym, Pethick, & 
Sutherland (1971) for the outer crust (p < 0.001 fm -3 ). 

Simulations of supernovae explosions (Burrows & Lattimer 
1986; Pons 1999) show that the protoneutron star has neither an 
isentropic nor an isothermal profile. For simplicity we will as- 
sume a constant temperature inside the star and attach a cold 
crust for the outer part. This schematizes the temperature pro- 
file of the protoneutron star. 

4. Results and discussion 

4.1. Composition of stellar matter 

First we will discuss the populations of beta-equilibrated stellar 
matter, by solving the chemical equilibrium conditions given 
by Eqs. dl5> . supplemented by electrical charge neutrality and 
baryon number conservation. In Fig. [5] we display the relative 
particle fractions as a function of the baryon density for several 



values of the temperature T — 0, 10, 30, and 50 MeV. The 
upper panels show the particle fractions when stellar matter 
contains only neutrons, protons, electrons, and muons, whereas 
the lower panels include the appearance of £~ and A hyperons 
treated as noninteracting particles. 

First let us discuss the purely nucleonic case. We notice that 
temperature effects influence the populations mainly in the low- 
density region. In fact, the presence of tails in the Fermi distri- 
bution makes it possible to create particles at any density and 
thus typical production thresholds, like for muon creation, dis- 
appear at finite temperature. On the other hand, thermal effects 
are less important at high density, where the particle popula- 
tions do not change appreciably when increasing temperature. 

The same considerations hold true when hyperons are in- 
cluded (lower panels). In this case, we notice that the YT and A 
thresholds disappear and that hyperons become more and more 
abundant at low density with increasing temperature. The lep- 
ton fractions at low density amount to about 10% of the total 
particle population if the temperature is high. We notice that 
the appearance of the IT induces a rapid deleptonization, no 
matter the value of the temperature. This happens because it is 
energetically more convenient to maintain the charge neutral- 
ity through IT formation than /3 decays. The stellar core turns 
out to be mainly populated by neutrons and protons, with an 
appreciable fraction (« 30%) of hyperons. 
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In Fig. |4] we show the particle fractions as a function of 
baryon density, for different values of the temperature, in the 
case that neutrinos are trapped. As in the previous figure, the 
upper panels show the particle fractions when stellar matter 
contains only nucleonic species, whereas the lower panels in- 
clude the appearance of noninteracting X~ and A. The elec- 
tron fraction is larger in neutrino-trapped matter than in the 
neutrino-free case; therefore the proton population is larger as 
well, and altogether the resulting EOS will be softer. The ap- 
pearance of muons is shifted to higher density values, because 
the onset is now determined by the difference between the elec- 
tron and the neutrino chemical potentials, fi e - yi Ve . However, if 
the temperature increases, muons appear at lower density. 

Neutrino trapping has also strong consequences for the on- 
set of hyperons. In fact, the onset of the Yr is shifted to high 
density, whereas A's appear at slightly smaller density. This is 
due to the fact that the E~ onset depends on the neutron and lep- 
ton chemical potentials, i.e., pt„ + fj. e -fi Ve , which stays at larger 
values in neutrino-trapped matter than in the neutrino-free case, 
thus delaying the appearance of the £~ to higher baryon density 
and limiting its population to a few percent. On the other hand, 
the A onset depends on the neutron chemical potential only, 
which keeps at lower values in the neutrino-trapped case. We 
notice that the appearance of the 2~ hinders muon production, 
which are not present at all. When the temperature increases, 
more and more hyperons are present also at low densities, but 
they still remain a tiny fraction of the total baryon density in 
this region of the protoneutron star. Altogether, the hyperon 
fractions are much smaller than in the neutrino-free matter. 



Fig-Elshows f or comparison the particle fractions with in- 
teracting hyperons in cold beta-stable matter with and without 
neutrinos. With respect to the previous results assuming nonin- 
teracting hyperons, one observes in neutrino-free matter a shift 
of the A threshold to substantially lower density, whereas the 
2T onset remains practically fixed. This can be related to the 
properties of the hyperon single-particle potentials at the rele- 
vant densities (Baldo et al. 1998, 2000a). In neutrino-rich mat- 
ter, the repulsive nature of the interactions at high density leads 
now even to a complete suppression of the S~, and only the A 
is present in the matter. 

4.2. Equation of state 

The resulting EOS is displayed in Fig. |6j where the pressure 
for beta-stable asymmetric matter, with and without neutrinos, 
is plotted as a function of the baryon density at temperatures 
T = 0, 30, and 50 MeV. Let us begin with discussing the case of 
neutrino-free matter, shown without (upper left plot) and with 
hyperons (lower left plot). We notice that thermal effects pro- 
duce a slightly stiffer EOS with respect to the cold case, and 
that at very high densities they almost play no role. The inclu- 
sion of hyperons produces a dramatic effect, because the EOS 
gets much softer (slightly less with interacting hyperons due to 
the repulsive effect of interactions at high density), no matter 
the value of the temperature. We observe, for the T = case, 
two wiggles associated with the onset of the E~ and A hyper- 
ons, which however tend to disappear when temperature is in- 
troduced, due to the disappearance of the thresholds, as shown 
in Fig. 13 

In the right-side panels we show the corresponding 
neutrino-trapped case. The EOS is slightly softer than in the 
neutrino-free case if only nucleons and leptons are present in 
the stellar matter. Again, the presence of hyperons introduces 
a strong softening of the EOS, but less than in the neutrino- 
free case, because now the hyperons appear later in the matter 
and their concentration is lower. A single large wiggle is now 
present at T = in the density range where both E~ and A ap- 
pear. Thermal effects are rather small also in this case, except 
for the disappearance of the hyperon onsets. 

4.3. Global stellar structure 

In Fig.^we show the gravitational mass (in units of the solar 
mass M Q ) as a function of the central baryon density for stars 
without (upper panels) and with (lower panels) hyperons. In the 
former case we notice that both finite temperature and neutrino 
trapping reduce slightly (by about 0.1 M ) the value of the 
maximum mass, which is about 1 .8 M for a cold neutron star. 
This is due to the fact that at high baryon density the EOS gets 
softer with increasing temperature. This behavior is at variance 
with the results reported in Prakash et al. (1997), where the 
critical mass increases with increasing entropy and, therefore, 
also with increasing temperature. 

In this regard, we should notice that in our calculations we 
are considering an isothermal profile, whereas in Prakash et 
al. (1997) the profile is isentropic. However, this cannot be the 
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Fig. 6. Pressure as a function of baryon density for beta-equilibrated matter at temperatures T — 0, 30, and 50 MeV. The left-hand 
(right-hand) panels show the EOS in neutrino-free (neutrino-trapped) matter, with nucleons only (upper panels) and nucleons plus 
hyperons (lower panels). The lower panels show results with interacting hyperons at zero temperature (faint lines) and with free 
hyperons at different temperatures. 



origin of the observed difference in the change of the maxi- 
mum mass as a function of temperature and entropy. In fact, 
the protoneutron star stability is mainly governed by the EOS 
in the central high-density part, and therefore an increase of 
entropy or temperature should be equivalent. Thus the value of 
the maximum mass at a given temperature in an isothermal cal- 
culation could be traced back to the corresponding isentropic 
calculation just from the entropy-temperature relationship for 
the EOS at the high density present in the core of the star. A 
more likely explanation can be found in the behavior of the 
EOS with increasing temperature in our microscopic calcula- 
tions. In fact, the decrease of maximum mass with temperature 
is in our case due to the decrease of asymmetry and to the con- 
sequent softening of the EOS. It has to be noticed that for a 
completely local interaction there is no dependence on tem- 
perature of the EOS except for the kinetic part. In Prakash et 
al. (1997) the interaction is indeed mostly local, with only a 
non-local correction, whose contribution is unfortunately not 
separately reported. In our calculations the whole interaction 
part is temperature-dependent, due to the intrinsic non-locality 
of the single-particle potentials, and therefore the temperature 
dependence can be quite different and larger. 



The lower panels of Fig.0show the configurations of stars 
containing hyperons. We observed in Fig.[6]that the EOS soft- 
ens considerably when hyperons are included, both in neutrino- 
free and neutrino-trapped matter. As a consequence the mass 
- central density relation is also significantly altered and the 
value of the critical mass decreases by a large amount, down 
to about 1.3 M for neutron stars and 1.5 M Q for protostars. 
Since the former value falls below the mass of the best observed 
pulsar, i.e., 1.44 M Q , the EOS of high-density nuclear matter 
comprising only baryons (nucleons and hyperons) is probably 
unrealistic (even taking into account the present uncertainty of 
hyperonic two-body and three-body forces) and must be sup- 
plemented by a transition to quark matter. This has been dis- 
cussed extensively in (Burgio et al. 2002a, 2002b; Baldo et al. 
2003; Maieron et al. 2004). 

Nevertheless, it is interesting to observe a behavior opposite 
to those of nucleonic stars, which could even survive the inclu- 
sion of quark matter, namely the maximum mass of hyperonic 
protostars is larger (by about 0.3 M Q ) than the one of cold neu- 
tron stars. The reason is the minor importance of hyperons in 
the neutrino-saturated matter, which leads to a stiffer equation 
of state, see Fig. |6] and to a larger maximum mass. This fea- 
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Fig. 7. The gravitational mass (in units of the solar mass M Q = 1.98 x 10 33 g) as a function of the central baryon density 
(normalized with respect to the saturation value po = 0.17 fm~ 3 ) at temperatures T = 0, 30, and 50 MeV. The left-hand (right- 
hand) plots regard neutrino-free (neutrino-trapped) matter. The upper plots describe configurations without hyperons, while in 
the lower plots the thin solid curves denote configurations of cold stars employing interacting hyperons, whereas the remaining 
curves show warm stars at different temperatures using noninteracting hyperons. 



Table 1. Characteristics of the maximum mass configurations 
for different stellar compositions and temperatures. 



Composition 


T (MeV) 


M/M Q 


R(km) 


Pc/P0 







1.86 


9.5 


8.2 


N,l 


10 


1.82 


9.5 


8.1 




30 


1.73 


9.7 


7.7 







1.76 


9.1 


8.8 


N,l,v 


10 


1.75 


9.2 


8.5 




30 


1.65 


9.5 


8.1 







0.93 


10.2 


7.8 




10 


0.95 


11.0 


6.7 




30 


1.00 


13.0 


4.8 


N,H mt J 





1.25 


8.8 


11.5 







1.34 


10.6 


6.9 


N,H h<x ,l,v 


10 


1.35 


11.0 


6.1 




30 


1.35 


12.2 


5.5 







1.48 


9.6 


8.8 



ture could lead to metastable stars suffering a delayed collapse 
while cooling down, as discussed in (Prakash et al. 1997; Pons 
et al. 1999). Our main results for the maximum mass configu- 
rations are summarized in Table 



5. Conclusions 

In this paper we have studied the structure of (proto)neutron 
stars on the basis of a microscopically derived EOS for bary- 
onic matter at finite temperature. Configurations with or with- 
out trapped neutrinos were considered. It was found that for nu- 
cleonic stars the thermal effects systematically reduce the max- 
imum masses with increasing temperature. This effect is argued 
to be related to the strong non-locality of the mean fields as 
obtained from the microscopic calculations. Neutrino trapping 
further reduces slightly the maximum mass, as a consequence 
of the additional softening of the EOS, but the dominant effect 
appears to be produced by the temperature, if values up to 30 
MeV are considered. These results indicate that the maximum 
mass of a nucleonic neutron star could be determined not by its 
mechanical stability but rather by the stability of the hot pro- 
toneutron star progenitor. 

On the contrary, we find that the maximum mass of a hy- 
peronic protostar is substantially larger than the one of the cold 
star, because both neutrino trapping and finite temperature tend 
to stiffen the EOS. Trapping shifts the onset of hyperons to 
considerably higher density and reduces their concentrations, 
in particular the E~ hyperon appears only in the higher density 
region or disappears completely. 
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Table 2. Parameters of the EOS fits, Eqs. (I25I26> . for symmet- 
ric nuclear matter (SNM) and pure neutron matter (PNM). 





0[ 


(22 


b 


h 


b 2 


CO 


Cl 


d 


E/A, SNM 


105 


74 


-473 


-464 




586 


381 


1.26 


El A, PNM 


109 


64 


34 


-240 




249 


164 


1.97 


F/A, SNM 


41 


116 


-180 




-174 


293 




1.57 


F/A, PNM 


21 


116 


101 




-131 


191 




2.62 



However, as in the case of cold neutron stars, the addition 
of hyperons demands for the inclusion of quark degrees of free- 
dom in order to obtain a maximum mass larger than the obser- 
vational lower limit. A consistent treatment of this aspect is left 
to a future work. 



Appendix A. Parametrizations of internal and free 
energy 

For the determination of the composition of beta-stable matter 
and for the solution of the TOV equations, it is very useful to 
provide analytical fits of the internal energy E /A(T,p, x p ) as 
well as the free energy F/A(T,p, x p ). 

It turns out that for both quantities the dependence on pro- 
ton fraction can be very well approximated by a quadratic de- 
pendence, as at zero temperature (Bombaci & Lombardo 1991; 
Baldo et al. 1998, 2000a): 

|(r,p,*)« |(7\p,* = 0.5) + (l -2x) 2 E sym (T,p), (22) 

where the symmetry energy E sym can be expressed in terms of 
the difference of the energy per particle between pure neutron 
(x = 0) and symmetric (x = 0.5) matter: 



1 d(E/A) 



4 dx 

* f(7>,0)-f(r,p,0.5). 



(23) 
(24) 



Therefore, it is only necessary to provide parametrizations of 
both quantities for symmetric nuclear matter and pure neutron 
matter. We find that the following functional forms provide ex- 
cellent parametrizations of the numerical results: 



(T,p) = (flit + a 2 t 2 ) + (b + feiOp + (c + ci)p d , (25) 
(T,p) = (a l t + a 2 t 2 )ln(p) + (b + b 2 t 2 )p + c l) p d , (26) 



where t = 77(100 MeV) and E, F and p are given in MeV and 
frrT 3 , respectively. The parameters of the different fits are given 
in Table|2l 
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